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XIX. On a New Method of Approximation applicable to Elliptic and Ultra-elliptic 
Functions. — Second Memoir*. By Chaeles W. Mebbifield. Communicated by 
W. Spottiswoode, Esq., F.B.S. 

Eeceived March 20, — Eead April 3, 1862. 

Since my first memoir on this subject was read before the Society, Mr. Sylvester has 
published a method, more general than mine, of applying rational approximation to 
facilitate the computation of the integrals of irrational functions. This method, at 
which he had arrived independently, included, a majori, the one which was the subject 
of my memoir. Aided by his papers, my subsequent studies have enabled me to view 
the method with more generality, as well as with more precision and completeness of 
detail, and I am now able to present it in a sufficiently finished and practical form for 
the immediate use of the computer. I have also computed auxiliary Tables, to render 
its application easier in certain cases. 

Any rational formula, which gives approximately the value of a function to be inte- 
grated, may be integrated in lieu of it, and the result will in general be an approximate 
value of the integral sought. But for such a process to be of any practical utility, the 
convergence of the formula must be excessive, for the complexity of the integral forms is 
so great that the labour would be enormous, unless the terms were very few in number. 
In the discovery of formulae sufficiently convergent for the purpose, lies the success of 
the method. 

We are by no means restricted to functions under a square root, or even to pure 
radical forms at all. The principle applies with equal generality to functions which are 
given implicitly as roots of equations, and thus to a class of differential equations ; and 
Mr. Sylvester has well remarked that these formula not only afford facilities for com- 
putation, as by a method of quadratures, but also enable us to assign superior and infe- 
rior limits to an integral, without losing its generality of form. 

I shall begin with the approximation to the square root, giving it in its general form, 
and explaining its exact analytical signification. I shall then show its application to 
Elliptic Functions, and how, in the ordinary cases, certain simple reductions can be 
effected, which greatly lessen the labour of computation ; and I shall give these reduc- 
tions for the cases more commonly occurring, with some examples and working formulae. 
I shall then add a short account of the extension of the method. 

The paragraphs in the first two sections of this paper bear a consecutive number for 
convenience of reference. 

* For the First Memoir, see the Philosophical Transactions for 1860, p. 223. 
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Section I. — Approximants to the Square Boot. 

1. Mr. Sylvester gives, for the approximants to the square root, the following state- 
ment ; — 

" Let r be an approximate value of \/N ; then by that mode of application of Newton's 
method of approximation to the equation a?=N, which is equivalent to the use of con- 
tinued fractions, we may easily establish the following theorem, viz., that 

r 2 + N r 3 + 3rN y 4 + 6r 2 N + N 2 r 5 + 10r*N + 5rN 2 

2r ' 3r 2 + N 5 . 4r 3 + 4rN ' 5r 4 + 10r 2 N + N 2 ' 

will be successive approximations to ^/N." 

2. Their general form is 

(r+VN)' + (r-VN) « N 
# — (r+VNf-tr-^N) 1 ^ iN > ........ (1.) 

which is always rational. In this form the approximation to y'N as i increases is 
obvious. The method of my previous memoir is simply the particular case of i=2 k . 

3. If we wish to approximate to N~*, we may take the reciprocal of (1.), or, what is 
simpler, we may divide (1.) by N, thus obtaining 

(r+yrN) l +(r--i/N) i 1 (9 , 

* = (r+VN)M^~VN)*7F * * * ^ 

Before we can integrate these formulae, we must reduce them by means of the method 
of rational fractions ; the simplest and most general way is as follows : — 

4. Let g be an ith root of unity ; then, obviously, 

log (1— a?)=log(l — g#)+log(l— fx)+ ..... +log(l-gfa). 
Multiplying the differential coefficient of this by (-—#), we obtain 



+ l~e 3 #+* * * * * -+i 



1— X % \—qx * \ — q l X l l~fx ' "••••• I 1— g*# 



and since ; .=1 + - . and -z =14.- ,, 



ii"i __ ,,2^-ri _^3^,+ • "f 



1— a?* 1-ga? ' 1 — g 2 a? j l — fx r ' ^1-gV 

Making #1=^^-^, we may thus divide N ± * into i fractions, each of the form 

I fr+4/N)+g*(r-*/N) 

f>+*/N)-g*(r--4/N) iN > 

A being any integer not exceeding i. 

5. If we add the pairs Jc and i—lc, we obtain for the sum of the pair, 

2 (r+y'N) 2 -2(r~^/N) 2 

(r-f */N) 2 + (r- y'N) 2 -^ + g*-*) (r 2 -N) ^ 

8rN 8r 

or 



2(r« + N)-(j*+g«-*)(r«-N) U1 same denominator' 
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according to whether the upper or lower sign be taken. Now, because g is an ith root 



2/c7T 



of unity, £*+£*-*= 2 cos — r-, and the sum of the pair reduces itself, for y'N, to 



y* 



i 

i 

1 
i 



4rN 



2rN 



(^ + N)-cos -^ (r 2 ~N) * r 2 sin 2 ^ + Ncos 2 -^ 

2 V y Z 2 



2rN 



1 



2rN 



i- 



N + sin 2 -^ (r 2 - N) * r 2 - cos 2 -* (r 2 - N) 



• (30 



/ 



For N k we have the simpler forms, 



Zh— — ~ 



4r 



1 



2r 



'fc 



* ( r 2 + N)~cos ?*? ( r s_ N) * ^sin^^ + Ncos 2 ^ 



1 

* 

I 






2r 



N + sin 2 ^(r 2 -~N) * r 2 



cos 2 — (r 2 — N) 



(4-) 



/ 



All that remains is to integrate these terms, and sum them. 

6. Our grouping the terms in pairs has limited the value of k to range from i to ^ ($— 1) 
when i is odd. There is an odd term which, however, presents no difficulty, being 

simply - in the case of v^N, and -^? in the case of N"~*. When i is even, k is limited to 



r 2 + N r 2 4- N 

range from 1 to \i— 1, and the odd term becomes — : — in the case of y'N, and —^r- 

in the case of N"*. It is important to bear in mind that the term just mentioned is an 
odd term, and therefore not affected with the coefficient 2, which appears in the terms 
composed of pairs corresponding to imaginary roots. 

7. The value of % which I consider to be most useful for general purposes, is i=8 : in 

this case the odd term becomes — z — or Q1VT , and the other values of — are three in 

number, viz. 22° 30', 45°, 67° 30'. With proper precautions i=8 will almost always 
give seven or more figures correct. 

8. If we now give infinite values to Jc and i and pass from the summation to the defi- 

/ . F 

nite integral, we have (putting A=t 



2 = N""*: 



i 



& 



2rd\ 



r 2 + (N— r 2 )cos 2 A7r 



__2 r C** d$ 

•Jo i + sm a^ 



and since 



C** d<p \v 

J H-^.sm 2 <j>~"" V\ + ~p 



this is an identical equation, as it ought to be. 

9. This use of approximants, therefore, is simply the application of the method of 

3 l2 
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quadratures to a definite integral, which we substitute for the surd proposed for 
evaluation. 

10. It would appear at first sight that a full application of the method of quadratures 
in the ordinary way, with the help of differences, would give better results than the 
mere summation of the ordinates. But this is not the case ; for the differences diverge 
immediately. If we use differential coefficients for the quadrature, instead of differences, 
we have an opposite anomaly, namely that the correction of the summation appears to 
be absolutely nil, inasmuch as the differential coefficients which appear in the series are 
all of odd order, and the numerator of each of them contains the factor sin <p cos <p, which 
vanishes at both the limits and -Jt. Legeto&e has discussed this point. See the 
Appendix to the second volume of his ' Fonctions Elliptiques,' p. 578. 

11. The application of the method to integrations, then, lies in the substitution for 



Jl 



* M nj „ C l f* 2Mr.dk. dt 

at of 



V^N u ° U1 J J ^+(N~r 2 )cos 2 A7r 5 

in which, since X and t are perfectly independent of each other, we may change the 
order of integration, thus obtaining 

'* f C l 2Mr . dt 



i 



d\ ; 



, o r 2 + (N~r 2 )cos 2 A7r 

and the rest of the operation depends upon our being able to perform the integration 
in } generally, and then to determine the integral in X by quadratures. The great 
advantage of the method turns upon the easy application of the method of quadratures, 
in consequence of our not requiring to difference the ordinates. 

12. One way of exhibiting generally the degree of convergence is as follows : 
N H always lies between 

(r+i/N)«+(r-*/N)« (r+ */N)*-fr- */N)< 

and the error of either is therefore always less than their difference, 

^ (r+^/N) 2i -(r->/N) 2i 

13. There is another mode, by which, in any given case, we may see how far it is 
necessary to carry our work in order to obtain a given number of decimals correctly in 
the result. Let 6 m be determined by the equation 

Jo cos« OT — m J cos 9/ 

/N r 

and let sin 1 = ^-7- or — t^t, whichever may be less than unity ; then the mth approximant 

will be ^~j~. This is easily seen from the general term of the approximant, since 



j 



J9 1 / 1-f sinfl ' 

cos0-~2 lo £« ll^-sinS 
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14. A table of meridional parts, such as is given in the books on Navigation, if carried 
far enough, would solve this equation. I have calculated an auxiliary Table for the 
purpose, as follows : — 



Let cosec <p— 1 =£, log f tan / - w + « <P ) =fr then 






,s z 



:^l0g2^^l0g^+ll0g(l+^) 



log2--logs+j-^+i£-i^+ 



2 



4 16 ' 48 128 



To bring this formula to the same unit as the common Table of meridional parts, 
we must multiply it by the number of minutes in the arc equal to unity, or by 
L=3437'74677 07849 4, whence we have \ L log, 2 = 1191-43224 08243 2, and 
iLlog. 10=3958-85223 39129 100. These data give the following Table, the argument 
being the common logarithm of z with its sign changed ; that is, the number of places 
which are correct : 



—\ogz. 


y- 


1 

2 
3 
4 
5 


5234-14859 

9117*70966 

13068-84816 

17026-92712 
20985-70200 



— log z. 


y- 


6 

7 

8 

9 
10 


24944-54650 
28903-39796 
32862-25012 
36821-10235 
40779-95458 



—log z. 


y- I 


11 

12 
13 
14 

15 

i 


44738-80681 
48697-65905 
52656-51128 
56615-36352 
60574-21575 



a/N 

15. As a simple example, let N= 3, r= 2; .*.- — = sin60°: the meridional parts for 

60°= 4527 ; and in order that the error may not exceed unity in the tenth place of 

40780 
figures, we must have m or i= ~~-- = 9 ; so that we must make i=9 at least, for the 

10th figure to be correct. 

16. These methods of course only exhibit the degree of approximation on the surd 
itself. The proportionate approximation is generally greater on the integral than on the 
simple surd, because the first approximant is usually so chosen as to be identical with 
the surd at one of the limits, and it is only near the other limit that the discrepancy 
tells. 

Section II. — Details of Reduction and Computation. 

17. The chief assistance, which can be provided a priori for the computer, consists in 

C t 2lS/Lrdt 

the exhibition and discussion, for the ordinary forms, of the integral \ ^rmZ^y^^ 

and of the auxiliary functions which present themselves in its reduction. 

18. In applying these methods to elliptic integrals, the radical and the first approxi- 
mant r must both be of a simple form, and it is advisable thatr 2 — N or N— r 2 should be 
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j 



of a square form. For the common form of the elliptic radical V^l— sin 2 6 . sin 2 <p) ? our 
choice is practically limited to 

(1) r=l, (2) r=cosip, (3) r=sin#. cosp, 

(4) r=cos#, (5) r=cos 6 . sin <p. 

And on these suppositions I now proceed to the integration of the general form of the 
reduced approximant for (J(l— sin 2 ^.sin 2 ^)~^<p=fe^<p. I omit mention of the con- 
stants of integration, because very slight changes in the function may alter them. The 
first of our three cases require, as they stand, no constant, and these are the most useful 



(^ctse^. 



(1) r=l, r 2 — N=sin 2 6. sin 2 9, 

2 



Zi 



k 



TT 



1 — sin 2 . cos 2 -r- . sin 2 

1 



)z h d<p = 2 / 1 — sin 2 6 . cos 2 -j J tan" 1 



1 — sin 2 6 . cos 2 — ) tan <b 



(2) 



T: 



Z 



cos?), r 2 — N=-~ cos 2 0. sin 2 CD, 

2 cos^ 



k 



2cos^> 



cos 2 <p -j- cos 2 . cos 2 — . sin 2 <p I—(l —cos 2 . cos 2 — - J . sin 2 <p 



{ 



$&d<P' 



( 1 — COS 2 . COS 2 — J loff. - 



#7r\* "1 



1 — cos 2 . cos 2 — j 



8 i / yforX^ 

! 1— sin<p( 1 — cos 2 0. cos 2 



;")'] 



(3) 



r: 



z 



k- 



sin . cos <p, r 2 — N= —cos 2 6, 

2 sin . cos <p 
( 1 — cos 2 . sin 2 -r- 1 — sin 2 . sin 2 <p 



J^<p = (1 — cos 2 6 . sin 2 



")" iog 1 (^ 



( 1— cos 2 .sin 2 -r-1 -f sin . sin <p 



2* ' «**Y • t • 
cos" 5 . sm J -7-I —sin . sin <p 



; 



(4) r=cos0, r 2 —- N=— sin 2 ^cos 2 <p, 



2 



2 



*■ 



1 + tan 2 . cos 2 — . cos 2 <t> 

1 T 



kv 



\jZ h dp=2 1 1 — sin 2 6 . sin 2 — ) tan~ J < 



cos#. tan <p( 1 



sin 2 0. sin 2 — 



'2 



(5) 



r: 



£fc : 



:cos0sinp, r 2 — -N=— cos 2 0, 

2 cos . sin <p 
cos 2 }■ . sin 2 © + cos 2 — . cos 2 
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ta k d$—2 1 cos 2 -y— cos 2 6 \ tan" 1 






'cos 6 / for . 

( cos 2 - — cos 2 6 

cos f y i 



COS 2 — COS 2 6 — ) l0£ . 



cos fl\ ^ / 

I eos<p/ s 2 #7r\i 

H 7 ( cos 2 $ — cos 2 — ) 

^ COS V\ I J 

(6) If we make £=tan -J<p, we obtain 

(1 - sin 2 6 . sin 2 p)-* ^=2(1 — 2 cos 2 . f+t)^ dt. 
Taking r=l— f, the terms which we have to integrate are of the form 

4(1 -t*)dt 

kit 
Putting ^ 2 =1 — cos 2 Q . cos 2 -r-, we have 

t 2 + 2qt+l 



y h d<p=- q log s (^ 



2qt+l 

The same expression serves for the integral 

2dt 



i/(l + 2cos2 0.^ + * 4 ) 



if we put q 2 = 1 — sin 2 # . cos 2 — . 

19. It will be observed that the first four cases, and the sixth, depend upon a radical 

of the form v"(l~ sin 2 A. sin 2 #), where oo is restricted to the selected values of -j. 

Assuming the modulus sin A not to vary, it would therefore in general be better to begin 
by computing the radical for the selected values. I have computed, and I append to 

this paper, a Table of this radical, the selected values of -j being 22° 30', 45°, and 67° 30', 

while A ranges by whole degrees from 1° to 90° inclusive. Every entry but the last in 
the 2nd, 3rd, and 4th columns of the Table was computed by myself in duplicate with 
Vega's ten-figure logarithms, by the help of two or more of the following formulae, some 
of which are from Legendre. 

20. Putting A for v"(l — sin 2 A. sin 2 ^), 

(1) Make sin A . sin a/=sin M ; then A=cos M, 

log sin M=log sin A+log sin a>, log A=log cos M ; or else 

(2) Make tan A . cos #=:tan M ; then A=cos A . sec M, 

log tan M=log tan A+log cos a>, log A=Iog cos A+ar. co. log cos M. 
Moreover, let L be the tabular angle nearest to the angtte M: it is not necessary to 
obtain the value of M : so that we have simultaneously, 

log sin M=log sin L+s, 
logtanM=logtanL±#, 
log cos M=log cos L+<? ; 
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then we shall also have, and with great approximation, 

log s=log (t . cos 2 L)+(£— # . cos 2 L) 

=log (c . cot 2 L) +(<?+<?• cot 2 L), 
log c=log (s . tan 2 L)+(s+s . tan 2 L) 

=log (t . sin 2 L)+(£— t . sin 2 L), 
log £=:log (s . sec 2 L)+(s—- s . sec 2 L) 

=log (c . cosec 2 L)+(c— c cosec 2 L). 

I have given the whole set of six, but my Table was computed, with the pair for log c. 
By way of example, I add a specimen copy of one of my working sheets. The use of so 
many as ten figures is not altogether unnecessary, because otherwise, when A is nearly 

equal to unity, the value of log (1— A) or of log frrK cann °t be had with exactness. 

21. The following formula will also be found in many cases preferable, both for 
exactness and facility, to the ordinary use of logarithmic tables by means of differences. 
These formulae, as well as those of the previous paragraph, are but applications of 
Taylor's theorem, reduced to a shape fit for the computer. Even where only seven 
figures are required their application is frequently much easier, and gives more exact 
results, than interpolation by differences. In what follows, x is supposed to be thf 
nearest tabular entry. 

22. To find log y from log tan y. — Let us assume simultaneously 

logy=logw+l, logtany=log tan#+#. 
Putting u=logx, £=logtan#, we have 

du sin 2x , d?u ^ T sin 2a? / n sin 2x 

Tz= ~W and d?= M ~2T V cos 2x ~ ~W 

M being the modulus of the logarithms. 
Hence, by Taylor's theorem, 



sin 2 t ( /Sin Qt \ i 

l=t — - — <1+Mn — g cos 2x) V nearly. 



Taking the logarithm, this becomes 

_ 7 , / , sin 2x\ , /sin 2x 

log Z=log(^ ^-J +t{-te—coa 2# 

, (t . sin x . cos x\ / . sin x . cos x , , rt , . 

=!og ( ~ x ) + -£— — ±t + 2t sm 2 x. 

The latter is the better shape for a working formula, because log sin x and log cos x 
are found in the same page and line as log tan x, while log sin 2x must be looked for 
elsewhere. The first term alone is sufficient when x is small ; but when x much exceeds 
45°, cos 2x changes its sign, and even the entire formula is insufficient. The maximum 
value of the coefficient of tin the second term is 1-0631, corresponding to #=78° 33' 26"*5. 
In many cases, where the first term alone is insufficient, a rough interpolation, made at 
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sight from the following Table, will answer the purpose ; it is a Table of the value of 
~o^r~— cos 2a; ) and of its logarithm, from #=45° to #=90°. 





45 


0*63662 


9-80387 


50 


0-73816 


9-86815 


55 


0-83149 


9-91986 


60 


0-91349 


9*96070 


65 


0-98041 


9*99141 



70 


1-02910 


0-01246 


75 


1-05658 


0-02390 


80 


1-06216 


0-02619 


85 


1-04334 


0-01843 


90 


1- 


0- 



The Table shows that, past 45°, the formula 



, , . j7sin#CQS#\ , 

log Z=log ( j +t 



X 



is a better approximation than when the +# is omitted. It is to be remarked that t is 
at its minimum for #=45°, and increases both towards #=0 and #=90°. Near the 
latter limit, where great accuracy is required, we must proceed as follows. 

Find the correction for the logarithm of the complement of the arc by the above 
process, and then find log (\ v—y) from log y. For this purpose, I observe that 
logy=log#+£ is equivalent to y=$.10 ± \ hence 

Now, let +A=10 ±! — 1, whence 

log (+mA)=log (+l)±^l— j2 M^ 2 nearly, and also 



i<>g(i»-y)=i«g(^-*)-(^F=F;-* M (?=j?; - • • • 

It is not often that the third term of either formula will be required. 

I have gone into all this detail, because the inverse tangent is continually presenting 
itself in all these integrations, and because no book that I know shows the proper way 
of handling it. 

23. The following constants are needed for these and similar formula? : — 



10+logm =9-63778 43113 00537, 
10+logl o =8-24187 73675 90828, 
10+logl' =6-46372 61172 07184, 
10+logl"=4-68557 48668 23541, 



logM=0-36221 56886 99463, 
■logl°=l-75812 26324 09172, 
log 1' =3-53627 38827 92816, 
log 1"= 5-31442 51331 76459. 



24. As an example of finding the inverse tangent, let it be required to find log^ and 
log (£*— y) from 

log tan#=9-02313 50437. Here we must take 
log tan #=9-02303 57359 

+t= 9 93078; .-. #=6° 1'10"=21670" 

~~~~ ~ fcr-aras 302330" 
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log sin#; 

log COS X: 
logt: 

ar. co. log#: 

, t sin 2os 

2nd correction 



9-0206346 
:9-9975988 
5-9969834 
0-9785662 

= 5-9937830 






log sin 3 #= 8-0413 

log 2=0-3010 

log £=5-9970 

log(2sin 2 a;)=4-3393 

2sin 2 #= 21,800 

— - — = 985,8 
2a? ' 



log 21670 
log 1" 



log a;: 






4-33585 89113 
4-68557 48668 
9-02143 37781 
+ 9.85786 



log #=9-02153 23567 



logZ=5-9937816 

+493 



17 



1007,6 

log mA= 5-9938309 2nd correction = 14,5 
log #=9*0214338 
ar. co. log (■gw— x) = 9*8339439 
log correction= 4*8492086 



The comma cuts off the 
eighth decimal. 



log 302330=5-48048 12441 

log 1 ;, = 4*68557 48668 

log (i7r-.r)=0*16605 61109 

correction = — 7 6 6 6 



log(-|7T-#)=0*16604 90443 



Verification. — The numbers corresponding to these logarithms of y and of \ir— y are 
0*10508 29743 and 1*46571 33525, the sum of which, to the very last figure, is 
exactly \w. 

25. To find log ~zr\ from logy. — Let log ^rj=log;^y+jp, and log ^==log #+£ ; 



y 

2qx \ 1 q 



y 



then log^=log [ ^ZTi )+2ff+~dri> nearly. This formula obviously fails where y is 



near unity ; in this case log -—? cannot be had with great accuracy, unless y itself be 

y 

given absolutely. All the cases of 1 2 _ -% may be included in the above formula by 

giving proper signs to p and q. It may save trouble to remark that # must not always 
be taken to the extreme limit of the Table, because log(#+l) and log (# — 1) have also 
to be taken out. As an example, let 



log 3/=0«-36290 63835 
log #= 0-36285 93030 
q- + 4 70805 
log 2=5-6728411 
ar.co.log(# 2 -l)= 9-3647540 

5-0375951 

log #=0-3628593 

log 2=0-3010300 

5-7014844 

2nd correction = —1225 



#=2-306, #+1=3-306, #-1=1-306 
log (#+l)=0-51930 28492 
log(#-l)=0-11594 31769 



sum=0-63524 60261 



2_ 



* 2 — 1 



i2— 



0-0001090 
135 



2nd corr n =0-0001225 



difFerence=0-40335 96723 
_ i?= _5 02761 

log ^4=0-40330 93962 



logp=5-7013619 
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This example has been so chosen as to admit of easy verification. In fact ^=2*30625, 

V + 1 529 

and log^—j==log 2^q=0 ; 40330 93959 24. The error is therefore only of three units 

in the tenth decimal place, where there was no reason to expect accuracy. 

26. The only other formulae which I shall give are the following, for finding the 
logarithm of a number, and vice versd. They are indispensable where more than seven 
figures are required. 

Let log (x+h)=log w+fc, then 

, T , /mh\ _1 mA , 

log £=iog ^— j +2 -£• nearl y> 

log h=log (Mwk)+ o k nearly. 

The values of log m and logM have been given in paragraph 23. 

27. As an example of the application of the method to the evaluation of elliptic 
integrals of the third class, let us take the integral 

^ dip 



X 



l (1 -sin 2 *. sin 2 ?) (l-sin 2 0.sin 2 <p)* 

for the values a=45°, 0=30°, p=60°. 

I have selected these values because they can be obtained without reduction or inter- 
polation from the Table of A(0, <p) which I have given, and also because sin 2 a = sin 6, 
and therefore the integral can be reduced to one of the first class, plus an inverse tangent, 
thus admitting of easy verification. For this case 



& k — — ^ --- / """ kit . \ 5 

(1 — sin 2 a.sin 2 <p) / 1 — sin 2 8. cos 2 -r-. sin 2 ? J 

2 sin k 1 

^z h d<p^z— — 1—. tan~* (cos cc. tan <p) 

J * Y • 2 • 2A 0™ C0Sa 



2sin 2 flcos 2 -r- 



sin 2 a — sin 2 fl . cos 2 — 

i 



(l~sin 2 0.cos 2 y ) tan""" 1 -! (l— sin 2 0.cos 2 y j tanp>. 



hit 
Making -j successively 22° 30', 45°, 67° 30', and, for the odd term, 90°, we find, after a 

few obvious reductions, that eight times the value of the integral is 



{ 



1 |+ An45 2o ,67|) + An45° ; 22i) }^5 5tan "^ COs45 °- tail60O} 



— 3^otan- l {cos 30°.tan 60°} -f A(30 o ?45 o ) tan-'{A(30°, 45°).tan 60°} 

_ tan-' { A(30°, g7j) .tan 60°}- .cos 2 22£ _ tan-' {A(30°, 22|) .tan 60°} .cos 2 67| 

A 2 (45° } 67J).A(30°,67|) A 2 (45°, 22i).A(30°,22|) 

3m2 
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As these inverse tangents range generally from 45° to 60°, I computed them by the 

(t sin ic cos#?\ 
— — 1 +t ; this being sufficient to 

give eight figures of decimals accurately. I found 

Iogtan- 1 -|cos45°tan60°£ =9-94747 15296, 
log tan- 1 {cos 30° tan 60°} =9-99246 23739, 

log tan- 1 -{A(30 o , 45°) tan 60°} =0-00766 92607, 
log tan-'{A(30°, 67|).tan 60°^ =9-99727 33807, 
log tan" 1 {A(30°, 221). tan 60°} =0-01665 09657. 

I hence obtained the following values : — 

Eor the positive terms. For the negative terms. 

7-10091 3039 1-13483 2441 

4-37212 6152 0-72539 4027 

2-70421 6251 1-66839 9478 



14-17725 5442 0-16728 4032 



3-69590 9978 3-69590 9978 



8)10-48134 5464 



1-31016 8183 value required 

A more exact value of the integral, otherwise obtained, is 

^F(30^60 o )+tan- 1 (~^-°j= 1-31016 8161, 

which differs from the previous value by 2 units in the eighth decimal place. 

28. In order to find how many places ought to have been accurately obtained, I 

13 

observe that the method followed gives N=yg, r =1, whence 

log (~ ^^=9-95491=% sin .64° 20' 30". 

The corresponding meridional parts are 5086*5, which must be multiplied by i=S, 
giving 40692*0. Eeferring to the Table in paragraph 14, I find that this nearly corre- 
sponds to ten places correct, and therefore that the integral ought to be correct to at 
least that extent. That it is not so, is due to my having curtailed the formula for 
finding the logarithms of the inverse tangents. But my object was only to give seven 
decimals correct, and my going beyond that was simply because, with a ten-figure 
Table, putting down the additional figures gave me less trouble (once I had to use more 
than seven) than abbreviation would have done. This remark may at first sight seem 
strange to any one who has not had some practice in using large Tables. But the loga- 
rithmic corrections are given in the shape of arithmetical complements; with reference 
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to the 10th figure, therefore, considered as an integer, the index is right as it stands, 
and we need not bestow thought on the proper placing of the correction, as we must if 
we use any other number of figures. 

29. If we had been content with five decimals, the calculation would have been very 
easy, for in that case we might have used six-figure logarithms, and have made i = 4, 
thus omitting the terms containing 22°J and 67°^. We should get 

7-10091 3039 1-13483 2441 

1-86022 6468 0-72539 4027 

4)5^24068 6571 1-86022 6468 

1-31017 1643 value required. 

30. It is worth while to notice a case which will sometimes occur, namely (using the 
notation of the last example), that the values may be so selected as to give, for one of 

the values of 1, sin a = sin 0. cos -j? and thus each of the terms into which §z k d<p was 

divided would become infinite. Of course the difficulty is only apparent ; for in this 

case the proper value is §z k dp = j/ 1 _ s i n g g f s i n « *)& °f which the integral may be at once 

found by differentiating the expression ^^ sin$ • 

Section III. — Extension of the Method. 

In respect of rapid approximation and precision of limit, the foregoing processes 
leave nothing to be desired, as far as concerns the radical of the square root ; but they 
do not go beyond that Mr. Sylvester has given an elegant extension of the method 
to radicals of a higher index, by means of symmetric functions * 

The more general problem before us is that of approximating to the integrals of 
irrational functions by means of rational substitutions. 

Let <p and $ be functional symbols, and y a function of z; then, that <p{z).y m and 
(p[z) : y m should both be approximations to <p{z\ depends upon y m approaching unity as m 
increases. Assuming that y m and y x are connected by the equation y m =^{m, y\\ our 
problem is to choose if so that, in the first place, the approximation shall be exceedingly 
rapid, and, in the next place, that <p{z).y m and (p(z):y m shall both (or at least one of 
them) be thoroughly manageable, and easily integrable. In the case of the approximants 

J fin s fill 

I am acquainted with three general methods which effect the object more or less. 
The first is the obvious one afforded by the Newtonian approximation to the roots of an 
equation ; viz., let a be a first approximate solution, obtained by trial, of the equation 
/a?=0, and call fx the differential coefficient of fx; then a second approximation is 

* See the Philosophical Magazine for December I860, Supplementary Numher, vol. xx. p. 525, note A. 



430 MB. C. W. MEKEiriELD ON A NEW METHOD OP APPROXIMATION 

a—jr^=b; a third approximation will evidently be b—^=c, and so forth. If we apply 

this method to the pure equation # w =p, the convergent terms which we obtain are as 
follows : — 

, (»— l)a n +p 
° na n ~ l ' 

__ m (n—l)\(n-l)a n +p\ n +n n .p.a< n - 1 ) „ 

— ri*a n ~ l {{n--l)a n +p\ n ~ l ' 

The second method is that of the reversion of series ; it is sufficiently discussed by 
Aebogast * 

The third method was suggested to me by Mr. Cayley's remark that Mr. Sylvester's 
third approximation is a particular case, for ft =2, of the common form (of the books on 

the binomial theorem) \/N= jz z jTjNf ■ L » i\ « • #> approximately, a being a first approxi- 
mation. In order to gain generality, and thereby symmetry, I shall pass from the par- 
ticular form \/N to the more general <p -1 N by the following Lemma: — 

and let x x , # 2 , # 3 , x 4 . . . . be determined by the system of equations, 

and so forth; also let N— a =|H>, and 

l—^,—^— ajf—a^— )(A +A 1 #+A 2 # 2 +A 3 # 3 + ), . . (3.) 

then a a t a 2 *„_! 

^1 — X ^ ^2—^9 ^3— A3 ^~ A w ' 

For, if we substitute these values in the equations (2.) after placing them in the 
following form, 

ju=:N— - %—a x X x -\ 

=a l oo 2 +a 2 x 2 x 1 I (4.) 

= #l#3 + #2#3#2 + #3#3#2#15 J 

and so forth, we obtain 

_ g i*o_gi*i +^Q_ fi? 1 A 2 4-a 2 A 1 + %A Q 
M ~ A x — A 2 — A 3 — CCC, (£).) 

which are the same equations as we should get by multiplying the two series in (3.) and 

equating to zero the coefficients of x and of its powers. The coefficient A =-, 
obviously. 

# Calcul des Derivations, pp. 288-296. 




(2.) 
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The coefficients A may now be found in a variety of ways ; by solving equations (4.) 
or (5.), by simple division, or by Akbogast's processes *. The object of the preceding 
lemma is to connect the quantities x n with the coefficients of division and of recurring 
series. Our results in any way are, 

a^ 4- «! V 



If for a 09 a„ &c. we substitute the coefficients of the binomial theorem, so as to make 



N =(«+#)*, we obtain 

a + a? »— (n+l)a" + (n-l)N * a ' 

Making %=2, we obtain Mr. Sylvester's approximants to the square root, and X w is then 
the coefficient of x n in the development by ascending powers of 

1 

(N — a 2 )— 2aa?— a?* ; ' 

and so far the method agrees with the Newtonian approximation by continued fractions ; 
but from this point the two methods diverge. For %=3, X n is the coefficient of x n in the 
development of 



and the successive approximants are 

2a 3 + N a 3 + 2N 4fl 6 +19fl 3 N + 4N 2 5a 9 + 45a 6 N + 30a 3 N 2 + N 3 

3a 3 * a > 2a 3 TN* a ' 10a 6 +16a 3 N+ W a > 15a 9 + 51a 6 N + 15a 3 N 2 ' a ' tea ; 

while the second approximant obtained by successive substitution is 

16a 9 + 51fl 6 N + 12a 6 N 2 +2N 3 

36a? + 36a 6 N + 9a 3 W ' a * 

What these methods all effect is simply a rational approximation to the value of y in 

du 

the equation <p(#, 3)=0. Then, making #==^, we have only to integrate in order to 
find the value of u. They thus constitute a means of approximately solving, in respect 
of u, differential equations of the form <p(z, j z ) = ; but they do not effect the solu- 

* See his ' Oalcul des Derivations/ pp. 26, 29 ; or De Mobgak, * Diff. Cale.' p. 331. 
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tion of this equation in respect of z, and still less do they solve the more general form 

It may suggest processes of reduction in some cases to remark, that there are many 
other functions of y m and <p(z)^ which will approximate to <p(z) as m increases, besides 
the simple product or quotient of <p{z) by y m . 

There is one point about these higher approximants, of which a solution, even if 
accompanied with considerable restrictions, would be extremely desirable, — I mean the 
resolution of the denominators into factors. I do not suppose that the problem, in its 
perfectly general form, admits of a compact solution ; but any class of cases, of even 
moderate generality, for which it could be elegantly solved, would probably have very 
useful applications. The criterion of convergence and the measure of approximation 
would also have their interest. 
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Specimen 


l sheet of work for the Table. 




Arc of 60°. 


log sin 22|= 9*58283 


96605 


83* 


log tan 2 L= 9-0913166 


log 


cos L= 


=9-97473 


27132 


log sin A = 9*93753 


06316 
02922 


96 
79 


log s= 


5-3160144 




c= 
log A= 


+ 


25545 


log sin M = 9*52037 


4-4073310 


=9-97473 


52677 


log sin L = 9-52039 

s = —2 


09944 
07021 




S + "I 

V -233 
s tan* L= J 










s tan 2 L = 


25546 




log c= 4-4073077 










log sin 45°= 9-84948 


50021 


68* 


log tan 2 L= 9-7781471 


log 


; cos L= 


=9-89794 


07883 


log sin A = 9-93753 


06316 


96 


logs= 4-1137763 




c= 


= — 


7797 


log sin M = 9-78701 


56338 
43344 
12995 


64 


3-8919234 

S + T 

} + 21 
s tan 2 L= J 




log A= 


=9-89794 


00086 


log sin L = 9-78701 
s = 4- 






s tan 2 L = 


7821 




log c= 3-8919255 










logsin67i= 9-96561 


53459 


21* 


log tan 2 L= 10-2501550 


log cos L= 


=9-77806 


24352 


log sin A = 9-93753 


06316 


96 


logs= 4-7816118 




c= 


= — 1 


07593 


log sin M = 9*90314 


59776 
99296 
60480 


17 


5-0317668 

s + i 

I +168 

s tan 2 L= J 




log A= 


=9-77805 


16759 


log sin L = 9*90313 

s = + 






s tan 2 L = 1 


07589 




log c= 5-0317836 










log tan 22^= 9*6 1722 


43146 


62* 


log sin 5 L= 8-6140815 


log cos #2| = 


=9-96561 


53459 21* 


log cos A = 9*69897 


00043 


36 


log t= 5-5451559 


log, 


sec L = 


= 0-00911 


84784 


log tan M = 9-31619 


43190 


08 


4-1591374 






9-97473 


38243 


log tan L = 9-31615 


92213 

50877 
14018 




t— i 


t 4-327 




c = 
log A= 


+ 


14430 


t = + 3 


tsin 2 L=J 

log c== 4-1591701 


= 9-97473 


52673 


tsin 2 L = 






log tan 45 =10-00000 


00000 


00* 


log sin 2 L= 9*3009948 


log< 


30s 45 = 


= 9-84948 


50021 68* 


log cos A = 






log t= 5-3429218 
4-6439166 


log 


sec L = 


= 0-04845 


06017 


log tan M = 9-69897 


00043 


36 


9-89793 


56039 


log tan L = 9-69894 


79790 
20253 
44047 




} +176 
tsin 2 L= J 

log c== 4-6439342 




c= 

log A= 


+ 


44049 


t = +2 


=9*89794 


00082 


tsin 9 L = 






log tan 67i= 10-38277 


56853 


38* 


log sin 2 L= 9-7730722 


log cos 67^= 


= 9-58283 


96605 83* 


log cos A = 9-69897 


00043 


36 


log t= 4-7719327 


log! 


sec L = 


= 0-19521 


55228 


log tan M =10-08174 


56896 


74 


4-5450049 






9-77805 


51834 


log tan L =10-08175 


16044 




I - 24 
t sin 2 L= J 

log c= 4-5450025 




c= 

logA= 


= — . 


35075 




59147 
35075 






tsin e L = 


=9-77805 


16759 



Mte.— The entries marked *, and the whole of the letter-press, were printed on the sheets. The letter A 
•n this page corresponds to 6 in the Table. 
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Table of the value of the function log A($, ») or log ./(l— -sin 2 0. sin 2 *) for four values 

of «, viz. 22° 30', 45°, 67° 30', and 90°. 



0. 


log A (*, 22° SO'). 


log A (*, 45°). 


log A (V, 67° 30'). 


log cos i. 


* 


o 

1 


9-99999 03138 


9-99996 69274 


9*99994 35386 


9-99993 38497 



1 


2 


9-99996 12643 


9-99986 77197 


9*99977 41349 


9-99973 53589 


2 


3 


9-99991 28792 


9-99970 24074 


9*99949 17311 


9*99940 44063 


3 


4 


9-99984 52048 


9*99947 10407 


9-99909 62308 


9-99894 07898 


4 


5 


9*99975 83051 


9-99917 36910 


9-99858 74990 


9*99834 42260 


5 


6 


9-99965 22633 


9*99881 04507 


9-99796 53618 


9-99761 43489 


6 


7 


9*99952 71805 


9-99838 14325 


9-99722 96070 


9*99675 07098 


7 


8 


9-99938 31764 


9-99788 67716 


9*99637 99831 


9-99575 27754 


8 


9 


9*99922 03891 


9*99732 66254 


9*99541 62004 


9*99461 99270 


9 


10 


9*99903 89748 


9*99670 11740 


9*99433 79300 


9*99335 14589 


10 


11 


9*99883 91084 


9*99601 06211 


9*99313 48042 


9-99194 65764 


11 


12 


9*99862 09825 


9*99525 51957 


9-99183 64168 


9*99040 43940 


12 


13 


9-99838 48090 


9*99443 51515 


9*99041 23213 


9*98872 39328 


13 


14 


9-99813 08170 


9*99355 07689 


9*98887 20337 


9-98690 41185 


14 


15 


9-99785 92542 


9*99260 23560 


9*98721 50300 


9-98494 37781 


15 


16 


9*99757 03868 


9*99159 02501 


9*98544 07479 


9*98284 16370 


16 


17 


9*99726 44986 


9*99051 48183 


9*98354 85857 


9*98059 63156 


17 


18 


9-99694 18915 


9*98937 64594 


9*98153 79027 


9*97820 63255 


18 


19 


9-99660 28857 


998817 56060 


9*97940 80200 


9*97567 00654 


19 


20 


9-99624 78189 


9-98691 27246 


9-97715 82198 


9*97298 58164 


20 


2\ 


9*99587 70469 


9*98558 83197 


9*97478 77462 


9-97015 17377 


21 


22 


9*99549 09429 


9*98420 29331 


9*97229 58056 


9*96716 58605 


22 . 


23 


9*99508 98979 


9'98275 71478 


9*96968 15661 


9*96402 60827 


23 


24 


9-99467 43204 


9*98125 15899 


9-96694 41603 


9*96073 01625 


24 


25 


9*99424 46358 


9*97968 69297 


9*96408 26837 


9-95727 57115 


25 


26 


9*99380 12870 


9*97806 38852 


9-96109 61968 


9-95366 01 869 


26 


27 


9-99334 47337 


9*97638 32245 


9*95798 37258 


9*94988 08840 


27 


28 


9-99287 54524 


9-97464 57677 


9*95474 42643 


9*94593 49269 


28 


29' 


9-99239 39363 


9*97285 23905 


9-95137 67741 


9-94181 92587 


29 


30 


9*99190 06948 


9*97100 40265 


9*94788 01877 


9*93753 06317 


30 


31 


9*99139 62526 


9-96910 16705 


9*94425 34100 


9-93306 55951 


31 


32 


9*99088 11517 


9*96714 63813 


9*94049 53211 


9-92842 04835 


32 


33 


9-99035 59484 


9-96513 92852 


9*93660 47788 


9-92359 14023 


33 


34 


9-98982 12144 


9-96308 15797 


9*93258 06231 


9-91857 42135 


34 


35 


9*98927 75363 


9-96097 45359 


9*92842 16784 


9*91336 45194 


35 


36 


9*98872 55150 


9*95881 95031 


9*92412 67607 


9*90795 76446 


36 


37 


9-98816 57654 


9*95661 79121 


9*91969 46797 


9*90234 86165 


37 


38 


9-98759 89157 


9-95437 12781 


9*91512 42488 


9*89653 21441 


38 


39 


9-98702 56072 


9*95208 12065 


9*91041 42888 


9*89050 25944 


39 


40 


9-98644 64943 


9*94974 93945 


9*90556 36388 


9*88425 39665 


40 


41 


9*98586 22425 


9-94737 76364 


9*90057 H640 


9*87777 98629 


41 


42 


9*98527 35294 


9*94496 78273 


9-89543 57674 


9*87107 34581 


42 


43 


9-98468 3 0429 


9*94252 19663 


9*89015 64024 


9*86412 74638 


43 


44 


9-98408 54812 


9*94004 21611 


9-88473 20869 


9*85693 40901 


44 


45 


9-98348 75524 


9*93753 06317 


9-87916 19193 


9*84948 50022 


45 
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Table 
(continued). 



J. 


log A (4 22° 30'). 


log A (*,45°). 


log A (0,67° 30'). 


log cos i. 


L 


d 


9-98288 79722 


9-93498 97136 


9-87344 50980 


9-84177 12731 


46 


47 


9*98228 74657 


9-93242 18620 


9-86758 09423 


9-83378 33303 


47 


48 


9-98168 67640 


9*92982 96543 


9*86156 89166 


9-82551 08951 


48 


49 


9-98108 66053 


9-92721 57944 


9-85540 86587 


9-81694 29168 


49 


50 


9-98048 77326 


9-92458 31150 


9*84910 00105 


9-80806 74967 


50 


51 


9-97989 08943 


9-92194 45794 


9-84264 30543 


9-79887 18039 


51 


52 


9-97929 68415 


9-91927 32846 


9-83603 81532 


9-78934 19787 


52 


53 


9-97870 63285 


9-91660 24627 


9-82928 59976 


9-77946 30249 


53 


54 


9-97812 01111 


9-91392 54820 


9-82238 76555 


9-76921 86852 


54 


55 


9-97753 89457 


9-9H24 58470 


9-81534 46333 


9-75859 13013 


55 


56 


9-97696 35878 


9-90856 71996 


9-80815 89399 


9-74756 16513 


56 


57 


9-97639 47918 


9-90589 33160 


9-80083 31625 


9-73610 87645 


57 


58 


9-97583 33090 


9-90322 81075 


9 79337 05490 


9-72420 97077 


58 


59 


9-97527 98869 


9-90057 56154 


9-78577 51006 


9-7H83 93361 


59 


60 


9-97473 52675 


9-89794 00084 


9-77805 16759 


9*69897 00043 


60 


61 


9-97420 01874 


9-89532 55788 


9*77020 61055 


9-68557 12291 


61 


62 


9-97367 53742 


9-89273 67330 


9-76224 53190 


9-67160 92909 


62 


63 


9-97316 15478 


9-89017 79885 


9-75417 74828 


9-65704 67649 


63 


64 


9-97265 94174 


9-88765 39622 


9-74601 21530 


9-64184 19615 


64 


65 


9-97216 96810 


9-88516 93618 


9*73776 04387 


9*62594 82593 


65 


66 


9-97169 30239 


9-88272 89739 


9-72943 51756 


9-60931 32999 


66 


67 


9-97123 01175 


9-88033 76506 


9-72105 11125 


9-59187 80116 


67 


68 


9-97078 16179 


9-87800 02961 


9-71262 51046 


9-57357 54170 


68 


69 


9-97034 81644 


9-87572 18497 


9*70417 63081 


9-55432 91617 


69 


70 


9-96993 03790 


9-87350 72689 


9-69572 63771 


9-53405 16846 


70 


71 


9-96952 88643 


9-87136 15099 


9-68729 96519 


9-51264 19176 


71 


72 


9-96914 42028 


9-86928 95088 


9-67892 33280 


9-48998 23640 


72 


73 


9-96877 69551 


9-86729 61579 


9-67062 76041 


9-46593 53400 


73 


74 


9-96842 76594 


9-86538 62846 


9-66244 57824 


9-44033 80750 


'74 


75 


9*96809 68303 


9-86356 46269 


9-65441 43168 


9-41299 62305 


75 


76 


9-96778 49569 


9-86183 58088 


9-64657 27859 


9-38367 51767 


76 


77 


9-96749 25025 


9-86020 43157 


9-63896 37732 


9-35208 80330 


77 


78 


9-96721 99032 


9-85867 44678 


9-63163 26324 


9-31787 89102 


78 


79 


9-96696 75670 


9-85725 03958 


9-62462 71226 


9-28059 88450 


79 


80 


9-96673 58731 


9-85593 60134 


9-61799 68925 


9-23967 02300 


80 


81 


9-96652 51705 


9-85473 49940 


9*61179 28070 


9-19433 24413 


81 


82 


9-96633 57778 


9-85365 07456 


9-60606 61106 


9-14355 53039 


82 


83 


9-96616 79819 


9-85268 63874 


9-60086 74357 


9-08589 44712 


83 


84 


9-96602 20377 


9-85184 47281 


9-59624 56795 


9-01923 45656 


84 


85 


9*96589 81672 


9-85112 82461 


9-59224 67793 


8-94029 60083 


85 


86 


9-96579 65594 


9-85053 90708 


9*58891 24439 


8-84358 45184 


86 


87 


9-96571 73697 


9-85007 89667 


9-58627 88989 


8-71880 01636 


87 


88 


9-96566 07187 


9*84974 93212 


9-58437 57166 


8-54281 91639 


88 


89 


9-96562 66934 


9-84955 11322 


9-58322 48116 


8-24185 53184 


89 


90 


9-96561 53459 


9-84948 50022 


9-58283 96696 


— log. infin. 


90 



